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Abstract: A growing body of evidence has been supporting the existence of so-called “dark molecular 
gas” (DMG), which is invisible in the most common tracer of molecular gas, i.e., CO rotational emission. 
DMG is believed to be the main gas component of the intermediate extinction region between Av^0.05- 
2, roughly corresponding to the self-shielding threshold of H 2 and 13 CO. To quantify DMG relative to 
HI and CO, we are pursuing three observational techniques, namely, Hi self-absorption, OH absorption, 
and TeraHz C + emission. In this paper, we focus on preliminary results from a CO and OH absorption 
survey of DMG candidates. Our analysis show that the OH excitation temperature is close to that of the 
Galactic continuum background and that OH is a good DMG tracer co-existing with molecular hydrogen 
in regions without CO. Through systematic “absorption mapping” by Square Kilometer Array (SKA) and 
ALMA, we will have unprecedented, comprehensive knowledge of the ISM components including DMG 
in terms of their temperature and density, which will impact our understanding of galaxy evolution and 
star formation profoundly. 
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1. Introduction 

The two relatively denser phases of the interstellar 
medium (ISM) are the atomic Cold Neutral Medium 
(CNM) traced by the Hi A21cm hyperfine structure line 
and the ‘standard’ molecular clouds (H 2 ) as traced by 
CO. CO is the most important tracer of molecular hy¬ 
drogen, which remains largely invisible due to lack of 
emission in the temperature range of molecular ISM. 
Empirically, CO intensities have been used as an indi¬ 
cator of the total molecular mass in the Milkyway and 
in galaxies through the so-called “X-factor” with nu¬ 
merous well-known caveats. Gases in these two phases 
dominate the masses of star forming clouds on a galactic 
scale. The measured ISM gas mass from Hi and CO is 
the foundation of many key quantities in understanding 
galaxy evolution and star formation, such as the star 
formation efficiency. 

A growing body of evidence, however, indicating the 
existence of gas traced by neither Hi nor CO. Com¬ 
parative studies (e.g. de Vries et al. 1987) of Infrared 
Astronomy Satellite (IRAS) dust images and Hi and 
CO gas maps revealed apparent ‘excess’ of dust emis¬ 
sion. The Planck collaboration (2011) clearly show ex¬ 
cess dust opacity (Fig.l) in the intermediate extinction 
range Av^ 0.05-2, roughly corresponding to the self¬ 
shielding threshold of H 2 and 13 CO. The missing gas, 
or rather, the undetected gas component is widely re¬ 
ferred to as dark gas, popularized as a common term by 
Grenier et al. (2005). They found more diffuse gamma- 
ray emission observed by Energetic Gamma Ray Exper- 
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iment Telescope (EGRET) than what can be explained 
by cosmic-ray H-nuclei interaction (H+X-factor*CO). 
Observation of the TeraHz fine structure C + line also 
help reveal dark gas in that the C + line strength in dif¬ 
fuse gas is stronger than what can be produced by colli- 
sional excitation with only Hi gas (Langer et al. 2010). 

A minority of the ISM communities argued that dark 
gas can be explained by underestimation of HI opacities 
(Fukui et al. 2014). Due to the limited scope of this pa¬ 
per, we will only discuss the dark molecular gas (DMG) 
hereafter, or more specifically CO-dark molecular gas. 

It is natural to infer from chemical and PDR models 
that molecular hydrogen would exist in regions where 
CO is not detectable. CO can be of low abundance 
due to photo-dissociation (Fig.2) in unshielded regions 
and/or can be heavily sub-thermal due to lack of colli¬ 
sion in diffuse gas. We strive to provide direct measure¬ 
ments and/or constraints of the physical conditions of 
DMG. Section 2 will focus on OH absorption. Section 3 
introduces a CO survey toward background continuum 
sources. Section 4 presents the combined analysis of CO 
and OH followed by a brief outlook of upcoming surveys 
in the final section. 

2. Where is the Hydroxyl? 

OH, or Hydroxyl, is the first interstellar molecule de¬ 
tected in radio bands (Weinreb et al. 1963). It can form 
quickly through a series of charge exchange reactions ini¬ 
tiated by cosmic ray once H 2 is present (van Dishoeck 
& Black 1988). OH can also form on grains. One of 
the main chemical paths associated with CO after OH 
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Figure 1. This figure is adopted from the Figure 6 in Planck 
(2011). Correlation plots between the gas column density as 
traced by [HI+Xco*CO] and dust optical depth at IRAS 100 
/im (upper left), HFI 857 GHz (upper right), 545 GHz (lower 
left) and 353 GHz (lower right). The color scale represents 
the density of sky pixels on a log scale. The blue dots show 
a binned average representation of the correlation. The red 
line shows the best linear correlation derived at low values. 
The vertical lines show the positions corresponding to Av = 
0.37 mag and Av = 2.5 mag. A single CO X-factor Xco = 
2.3 xlO 20 H 2 cm _2 /(K km s -1 ) was used. 

formation is 


OH + C + -> CO+ + H , 

(1) 

CO+ + H 2 -► HCO+ + H , 

(2) 

HCO+ + e“ -»■ CO + H . 

(3) 


We should expect wide-spread and abundant OH 
along with HCO + and C + . HCO + is accessible in mil¬ 
limeter bands. The main transition from C + is its fine 
structure transition in 2 THz band, which is impossi¬ 
ble to map from the ground. It is somewhat puzzling 
why large scale OH surveys of ISM has not been avail¬ 
able in the half a century since its discovery. In fact, 
thousands of hours of Arecibo time have been spent on 
searches for OH in galaxies with mostly negative results 
(e.g. Schmelz & Baan 1988). 

In contrast, Dickey et al. (1981) found OH in 
absorption against high Galactic latitude continuum 
sources. Important and extensive confirmatory absorp¬ 
tion measurements by Liszt & Lucas(1996) and Lucas 
& Liszt(1996) found that OH and HCO + are commonly 
observed against such sources. Lucas & Liszt(1996) 
found that ^30% of continuum sources having Hi in 
absorption exhibit HCO + in absorption. In light these 
results, the dearth of OH emission should be attributed 
to the excitation condition of OH rather than its abun¬ 
dance. 

The observed antenna temperature T A is 

T A = (T ex -T c )( l-e~ T ) (4) 

where Tq ~ 3.5 K is the continuum background temper¬ 
ature at L-band, comprising of the CMB and the galac¬ 
tic synchrotron emission. When T ex approaches To, 
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Figure 2. A schematic view of photo-dissociation region (Tie- 
lens 2005). It shows the locations of different transition lay¬ 
ers. We add the blue pane to indicate the possible location 
of OH. 


the apparent signal from certain line emission vanishes. 
Such gas, however, is suitable for absorption study when 
the telescope is trained toward background sources with 
Tc T ex as is the case when observing quasars and/or 
HII regions. 

Heiles & Troland (2003) published the well-cited Mil¬ 
lennium survey of 21-cm line absorption toward 79 con¬ 
tinuum sources. The ON-OFF technique and Gaussian 
decomposition analysis allow them to provide credible 
measurements of the excitation temperature and den¬ 
sity of Hi components spreading through the Milkyway. 
Unpublished OH absorption data were taken simulta¬ 
neously with the Millennium survey. Our preliminary 
analysis of these OH absorption data confirms the sus¬ 
picion that T ex of OH aggregates around the background 
temperature, and thus renders OH undetectable in emis¬ 
sion. 


3 . A Multi-transition CO Survey of 
Millennium Sources 

We conducted a follow-up CO survey of all Millennium 
sight-lines with OH absorption. Toward 79 published 
Millennium survey sources, 43 sight-lines exhibit OH 
absorption. These 43 sources with OH absorption have 
been observed in 12 CO J=l-0 and 2-1, 13 CO J=l-0 and 
18o J=l-0. The 8 sources with strong 12 CO J=l-0 and 
2-1 lines are also observed in the 12 CO J=3-2 transition. 

The J=l-0 transition of 12 CO, 13 CO and 18o were ob¬ 
served in March and May of 2013 and May of 2014 with 
the Purple Mountain Observatory Delingha (PMODLH) 
13.7 m telescope, Chinese Academy of Sciences. The 
spectra was obtained with position switch mode, the 
reference position is selected from IRAS Sky Survey At- 
la^] With 1000 MHZ bandwidth spectroscopy, the fre¬ 
quency resolution is 61 kHz, results in an approximate 
0.18 kms -1 channel width. 

The 12 CO J=2-l and J=3-2 data were taken with 
the Caltech Submillimeter Observatory (CSO) 10.4m on 
top of Mauna Kea in July, October and December of 

1 http : //irsa.ipac.caltech.edu/data/ISSA/ 
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Figure 3. The location of sources in galactic coordinate. Tri¬ 
angle represents the source without any absorption compo¬ 
nent in OH. Square represents the source with absorption 
component in OH. Yellow dot represents the source with¬ 
out any detection of CO transitions. Red dot represents the 
source with CO transition components, and these CO com¬ 
ponents correspond to all the absorption components in OH. 
Black dot represents the source with CO transition compo¬ 
nents, and these CO components correspond to some of the 
absorption components in OH, but there are no detectable 
CO transition corresponding to the rest of absorption com¬ 
ponents in OH. We call these kind of sources “part CO de¬ 
tection” . 

2013. The velocity resolution for 12 CO(J=2-l) spectra 
is 0.16 kms -1 . The velocity resolution for 12 CO(J=3-2) 
spectra is 0.11 kms -1 or 0.56 kms -1 due to a problem of 
spectroscopy. 

The distribution of these spectra in the galactic coor¬ 
dinate is presented in Fig. [3] 

The typical RMS of 12 CO J=l-0 observation is about 
0.06K, which corresponds to a CO detection limit of 
2.6xl0 13 cm -2 . 

The astronomical software package Gildas/CLASSrl 
was used for baseline removing, spectra combining and 
gaussian fitting. 

4. Comparison of Hi, OH, and CO 

We compare the Gaussian components seen in Hi ab¬ 
sorption, OH absorption, and CO emission. A total of 
115 Gaussian components were detected as specified in 
Heiles & Troland (2003). 52 such gas components have 
OH absorption. The majority of these 52 have CO emis¬ 
sion, except for 13 components, which are DMG candi¬ 
dates. There are no component with only CO emis¬ 
sion and no OH absorption. Three representative sets 
of spectra are shown in Fig. [4j Toward 3C192, only Hi 
is present, typical of CNM. Toward 3C154, there is a 
component with Hi, OH, and several CO and CO iso- 
topologue transitions, which should be representative of 
‘normal’ molecular clouds. Toward 3C409, there exists 
a component with Hi and OH absorption, but no CO 
emission in any transition observed. We posit that this 
is typical of DMG. The percentage of these three cate¬ 
gories are 55% CNM, 34% molecular clouds, and 11.3% 
DMG. 

2 http : / /www.iram.fr / IRAMFR/GILDAS/ 


Figure 4. Representative spectra. 3C192 sightline has only 
Hi seen in absorption. One component of 3C154 sightline 
has OH and Hi in absorption and CO and its isotopologues 
in emission. 3C409 sightline has one gas component with 
both OH and Hi , but no detectable CO transitions. 

The column densities of OH components were calcu¬ 
lated as 


Nqh 


8nkT ex uf m7 16 f 

—2 - 3j—~r / T 1667 dv 

Ai667 5 J 


(5) 


where Ai 667 = 7.778 x 10 _11 s _1 is the A-coefficient and 
T ex is its excitation temperature calculated based on a 
recipe similar to that for Hi absorption components in 
Heiles & Troland (2003). 

The CO column densities were calculated in two cate¬ 
gories. If only the J=l-0 transition of 12 CO is detected, 
the optical depth is assumed to be small and the exci¬ 
tation temperature is assumed to be the same as that 
of OH. If both 12 CO and 13 CO are detected, we derive 
the optical depth and the excitation temperature based 
on multiple transitions and Local Thermodynamic Equi¬ 
librium (LTE) assumptions. The recipe for deriving CO 
column densities can be found in Li (2002). 

The statistics of gas column densities (Fig. [5]) is con¬ 
sistent with the schematic picture presented in Fig. [3] 
and section 1. There is an apparent gas column density 
threshold for OH detection at around Av^0.05, above 
which OH and CO have similar distribution. OH turns 
out be a good tracer of diffuse gas with ‘intermediate’ 
extinction, namely, between the self-shielding threshold 
for H 2 and 13 CO. 


5. Discussion 

The expected location, abundance, and optical depth of 
OH should make it an excellent tracer of DMG. Due to 
insufficient collision in diffuse gas, however, OH is hard 
to detect in emission. This is likely the main reason 
why a galactic scale or even any large-scale OH map 
has not been accomplished. The realize its potential in 
quantifying dark gas throughout the ISM, the upcoming 
radio telescopes will be needed to conduct comprehen¬ 
sive absorption surveys. The Five-hundred-meter Aper¬ 
ture Spherical radio Telescope (FAST) is expected to 
start operation in late 2016. The unprecedented sensi¬ 
tivity of FAST and its early science instruments (Li et al. 
2013) should make feasible a Hi+OH absorption survey, 
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Figure 5. The histogram of Hi column density for the cold 
neutron medium(the blue filled rectangle), molecular gas(the 
green line filled rectangle) and dark gas (the red line filled 
rectangle). 

in the mode of the Millennium survey, but with 10 times 
more sources. The SKA1 will have the survey speed 
and sensitivity to measure gas absorption with a source 
density between a few to a few tens per square degree 
(McClure-Griffiths et al. 2014), which means that an all 
sky “absorption-image” is feasible and we will have ISM 
temperature and density everywhere! Based on simi¬ 
lar excitation and sensitivity considerations, ALMA is a 
powerful instrument to obtain systematic and sensitive 
absorption measurements of millimeter lines in diffuse 
gas. CO and HCO + in diffuse gas, in particular, will 
be much better constrained in terms of excitation tem¬ 
perature and column densities through ALMA absorp¬ 
tion observation than emission measurements. Com¬ 
bining both radio and millimeter absorption surveys in 
the coming decade, we will quantify DMG and provide 
definitive answers to questions like the global star for¬ 
mation efficiency. 
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